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A semiclassical Foldy-Wouthuysen transformation of the Dirac equa- 
tion is used to obtain the radiationless Bloch equation for the polari- 
sation density. 

1 Introduction and Motivation 

A complete picture of spin polarisation in accelerators and storage rings, either with 
or without synchrotron radiation, can only be obtained on the basis of evolution 
equations for combined spin-orbit distributions. See Article I. Moreover, if we are 
concerned with the effects of radiation, its simulation by classical white noise and 
damping does not suffice for all situations. For example we cannot obtain the Sokolov- 
Ternov effect by that means. In fact to include all the subtleties of radiation, a 
quantum mechanical approach is needed and then obtaining the 'complete picture' 
implies that we must begin by finding the equation of motion for the spin-orbit 
density operator in the presence of radiation. To ensure some level of transparency 
and trackability one begins by ignoring direct and indirect inter-particle effects so 
that at the classical level the beam would be described by a single-particle density 
depending on the six orbital phase space variables, the spin variables and on time as 
in statistical mechanics in '/i-space'. 

In the single-particle approximation, only positive energy two-component spin- 
orbit wave functions are needed. The appropriate quantum Hamiltonian is provided 
by a Foldy — Wouthuysen (FW) transform |IJ of the Dirac Hamiltonian and in 
order to get explicit results for time- dependent electromagnetic fields one has to use 
perturbation theory. Since we are interested in high energy behaviour in storage rings 
we do semiclassical perturbation theory, where the expansion parameter is Planck's 
constant, not 1/m etc. Before launching into the full blown calculation of the effects of 
radiation one should first obtain the transformed Hamiltonian for motion due to the 
Lorentz forces in the fields of the storage ring and then the corresponding equations of 
motion for the spin and orbital parts of the density operator. The evolution equations 
for the resulting classical distributions should then be derived. These are the tasks 
of this paper. Radiation will be considered elsewhere. 

a Updated version of a talk presented at the 15th ICFA Advanced Beam Dynamics Workshop: "Quan- 
tum Aspects of Beam Physics", Monterey, California, U.S.A., January 1998. Also in DESY Report 
98-096, September 1998. 
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The required Hamiltonian was already stated at first order in h by Derbenev and 
Kondratenko [0] as the starting point for their radiation calculations, but no construc- 
tion was given. In this paper we show how to obtain the transformed Hamiltonian to 
second order in H. Owing to space limitations we will be very brief but full details 
can be found elsewhere [[|; see also |4], ||, ^, [/J. 

2 The Dirac equation 

The underlying Dirac-equation is ih(dip/dt) = H op if) where: 



:= c(a op ■ n op ) + mc 2 ^ opA + e 
e(g-2)h 



op 

(lop,ii^op ■ B op ) + (7 op • E op )) . 



Pauli term 

We use the matrix representation of |§. 
3 The perturbation theory 

A unitary operator U op which performs the FW transformation (a 'unitary FW opera- 
tor') transforms a wave function ip in the Dirac representation into the corresponding 
wave function U op ip in the FW representation. In the FW representation the Dirac 
equation reads as ifo^- = H op ip where 

dU 

H op := U op H op U^ 1 + ih-^ U- 1 (1) 

and the electron wave functions, i.e. the 'positive energy wave functions' have only 
upper components which we denote by Xi;X2- Then in the Dirac representation the 
positive energy wave functions read as ip = U~p(xi, X2, 0, 0). Thus in this represen- 
ts 



tation the 'positive energy projection operator' A^ 5 reads as 



■■= U^l + lopA )U op . (2) 

By Eq. (^|) the FW transformation (and the FW Hamiltonian) is not unique because 
U op can be replaced by any unitary operator V op , such that t^opKp 1 i s even ('essential 
uniqueness' of the FW transformation). Our task is now to develop a perturbation 
theory to construct a U op . Then at k-th order (k > 0) in H, U op is approximated by 
an operator U op ,k- Starting with zeroth order, to satisfy Eq. (§), we must choose U op fi 
such that 

^ = Uo- p \ l -{l + l opA )U opfl . (3) 
Because for semiclassical perturbation theory we have||: 

AoV = - ( 1 + — (a op ■ nop) + -j-JopA) , (4) 

^ ^ J op J op ' 
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we can choose ||: 

U opfi := expf-i7 0Pi 47 p i5 arctan( — (S op • vr op ))) • 
In k-th order Eq. ([!]) reads as 

H op ,k '■= U ov hH ov U \ + ih — 7rr~U 



. dU p t k jj—i 



where -ff p,fc denotes the fc-th order approximation of the FW Hamiltonian. Having 
fixed Uop,o the remaining task is to choose the U op ^ for k > 1 such that H op ^ is even 
in fc-th order. That this procedure leads to a unique energy projection operator can 
be shown || by using a method due to Pursey |J. For 
k > 1 we compute the Z7pp,& recursively in the form 
U op ^k '■= U op ^U oPi k-i • • • Uop t iU oPj Q : where: 

U op ,j ■= exp(i[-^7 oPi4 , C oPi j-i]) , (j > 1) , 

with J op := a/^ 2 c 4 + c 2 (tt op ■ n op ) and where O op j is the odd part of the operator 
H op ,j, so that O oPt j anticommutes with 7 oPi4 and H op j — O op j commutes with 7op,4- 
Thus Oopj is the odd part obtained after the (j + l)-th step. 



In first order one gets H op i — || where: 

ijdk 1 e j. e ^ f mc2 9~^\ c$ B \ 

H op '■= ifiopAJop + ^(po P - -^X^j — I — ^—nopA^op ■ Hop) 



i" _ , „\ Top, 4^ ""op ' ^op)\/-'op ' ^ op) 



ehc(g — 2) 

8m J op (Jop + m c 2 
eH/g-2 1 mc 2 \ / - /w - A 

In second order, i.e. in H op2 , terms quadratic in E op and B op and gradient terms in 
E op and B op (e.g. the 'Darwin term') appear but no spin terms. 

4 The two-component formalism 

The two-component Hamiltonian ft,* the 'DK-Hamiltonian', is obtained from if* 
by the replacements 7 oPj4 — > 1 and S op — >■ a op so that: 

:= h d Q k p ,or b + \ ■ (^op • fiop) , (6) 

where h d ^ orb := J op + e<f) op and: 

- ._ e mc 2 p-2 ^ ec(g-2) 1 - 

op ' 2mc l J op + 2 j ° p + 4m J op ( J op + mc 2 ) 1 ° p ' ° pJ ° p 

+2^ (2^ + WpT^y J (7f - A Eop) + hx - • (7) 
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5 Nonrelativistic approximation 

Expanding fo* w.r.t. 1/m one gets in first order in 1/m ('nonrelativistic limit'): 
mc 2 + (l/2m)(7T op • 7T op ) + ecf) op — (ehg /4mc)(a op ■ B op ) which for g = 2 is the Pauli- 
Schroedinger Hamiltonian. 

6 The density operator in the two-component formalism 

The density operator p op reads as: 

1 

Po P - ^ 



i.Pop,orb @op ' ^iop) ) 



where p 0Pt0r b and £ op do not contain spin degrees of freedom. The normalisation of the 
density operator reads as 1 = Tr[p op ] = Tr[p 0Pt0r f,}. Expanding the Hamiltonian and 
the density operator to second order, the evolution equation for the density operator, 
the 'von-Neumann equation', reads in first order as: 



±_dpop l_r h d k] 

dt h [pop) opl 



Note that the second order parts of the Hamiltonian drop out of this von-Neumann 
equation because they are independent of spin. By Eq. (8) the von-Neumann equation 
is in first order equivalent to: 

= 7^ h y[Pop,orb > h^pfirb] + 2 ' ^°p) ~~ (^°P ' W)) ' (^) 

= + \&o P A 4) - \(C P A tiop) + \\U , hf pt0rb ] 

% ~* 
+ 2^-P°P' orb ' ^ op ] ' 

The terms proportional to Q op in Eq. ([|) account for the effect of the SG force on the 
orbital motion and the second and third terms on the rhs of Eq. (|T0D have the same 
structure as the Thomas-BMT equation. 

The terms h d ^ orh and Q op in the von-Neumann equation are not unique because 
the FW Hamiltonian depends on the chosen FW transformation. However, essential 
uniqueness allows the forms in Eqs. (|6|) and (|7|). 

7 The Weyl transform in the two-component formalism 

The 'Weyl transform' allows q-numbers to be represented by c-numbers. In the two- 
component formalism an operator K op is represented by its Weyl transform K wt via: 

K wttUX (r,P] t) := Tr[K op M op!V xA op (r,p)] , (v, X = 1, 2) , (11) 

where the operator A op is defined by: 

A op (r,p) := J d 3 x d 3 wexpQ[f- {r - r op ) +u- {p-Pop)]j , 



and where the operators -M oPi ii, ... are defined by: 

{M op ,vX x)» = <Wa , (a*, v, A = 1, 2) . 

Conversely one has: 

1 2 f 

K op = Yl M o P ,uX / d 3 r d 3 p K wt ,„xK P • (12) 

u,X=l 

Thus the Weyl transform is a 2 x 2 matrix valued phase space function; the t- 
dependence in K wt only occurs for time-dependent operators. In terms of its Weyl 
transform the trace of an operator K op reads as: 

Tr[K op ] = ^^tr[J d 3 r d 3 p K wt ) , (13) 

where tr denotes the matrix trace. 



8 The Wigner function 

The chosen normalisation of the density operator p op leads to: 



^tr[J d 3 rd 3 pp wt ], 



(14) 



for its Weyl transform p wt (see Eq. (0)) and one calls (l/8n 3 h 3 )p w t the 'Wigner 
function' of that state. In terms of Weyl transforms the expectation value of an 
operator K op reads as: 



< K o p >= 

87T 



^tr[ J d 3 r d 3 p p wt K wt ] . (15) 



9 The Wigner-Kirkwood expansion 

Since we are dealing with a beam, which in a high energy accelerator occupies a phase 
space volume many orders of magnitude greater then h 3 , we are very far from dealing 
with a pure state. Then, applying semi-classical perturbation theory to the density 
operator, its Weyl transform has the Wigner-Kirkwood' form: 

Pwt= s Po^ + h- Pi + fi 2 ■ p2 + (16) 

Classical part Quantum corrections 

where p , pi, p 2 , ... are of zeroth order in hf\. In reality, 'classical distributions' p do 
not exist so that one has to deal with Eq. ([16[) . 

b Note that p op has a related expansion and that this was used in Eqs. (9) and (10). 
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10 The Weyl transform of the Hamiltonian 

The Weyl transform of the Hamiltonian is hf^ t — /i* 6 + |(cr ■ f2 dfc ) where: 

hf rb :=J + e<P, 

=5 e mc 2 g ce(g - 2) 1 - 

^dfc := (— p + — -B + = 777- 2T ( ' B ) n 

mc J 2 2m J(J+mc z ) 

e (q — 21 mc 2 \ , w -=>. 



m V 2 J J( J + mc 2 



and where 7? := p — ~ c A and J := ^Jm 2 ^ + c 2 (tt ■ 11). 

11 The Weyl transform of the von-Neumann equation 

The Weyl transform of the von-Neumann equation is an evolution equation for p wt . 
In particular from Eqs. (Q) and fllOD and in first order using section 7, one gets: 

— P w *'° rb _|_ {h^ b , p wt ,orb} + ~{^dk,j , £wt,j} , (17) 

V v ' V v ' 

classical Liouville terms SG term 

= + {^orb 5 £wt} + &<dk A £^ + —{Qdk , Pwt,orb} , (18) 

V v ' S v ' 

Thomas— BMT terms SG backreaction term 

where { , } is the usual Poisson bracket w.r.t. r and p and where repeated indices 
are summed over. As in Eqs. (H) and fllO|) the second order parts of the Hamiltonian 
drop out. The vector £ wt is the polarisation density and Eq. (18) is its Bloch 
equation in first order ||. Note that these equations are not restricted to dipole 
and quadrupole magnetic fields. It is the p wt and £ wt which serve as the 'classical 
distribution functions' which we have been seeking. 

These equations are easily transformed from cartesian coordinates to 'machine 
coordinates' since with the Weyl transform one only has to deal with c-numbers 
instead of q-numbers |K| . After transforming to machine coordinates the zeroth order 
limits of the transformed Eqs. fllTD and ( |18| ) correspond to the classical Eqs. (17) and 
(40) in Article I. 

For FW transformations of the Dirac equation where time has been replaced by 



the longitudinal coordinate in a paraxial approximation see [11 



12 Radiation 



Now that the radiationless case is on a firm basis, one can include radiation. In a 
classical treatment of radiation effects one gets the Fokker-Planck and Bloch equations 
of Eqs. (24) and (39) in Article I. To include the Sokolov-Ternov effect one needs a 



full quantum treatment 12 
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